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THE FREQUENCY OF TORSIONAL VIBRATION OP A TAPERED BEAM 

By Robert P. Coleman 

SUMMARY 

A solution for the equation of torsional vibration o 
tapered beams has been found in terms of Bessel functions 
for beams satisfying the following conditions: (a) The 
cross sections along the span are similar in shape; and 
(b) the torsional stiffness of a section can be expressed 
as a power of a linear function of distance along the span. 
The method of applying the analysis to actual cases has 
been described. Charts are given from which numerical 
values can be immediately obtained for most cases of prac- 
tical importance. The theoretical values of the frequency 
ratio have been experimentally checked on five beams hav- 
ing different amounts of taper. 

INTRODUCTION 



The frequency of torsional vibration of a uniform 
beam can be calculated for a number of shapes of cross 
section. (See references 1 and 2.) The corresponding 
cases of tapered beams either have not been solved or are 
not readily available. A special case of longitudinal 
vibration of tapered beams treated by Nabl (reference 3) 
is applicable also to torsion for a special case (n a 2). 
In the present paper, the analysis has beon extended to 
the case of tapered "beans subject to rather broad condi- 
tions, and a closed solution of the differential equation 
of motion has been obtained in terms of Bessel functions. 
The result is expressed in such a form that the effect of 
taper is given independently of the effect of cross sec- 
tion. Hence, except for certain limitations imposed on 
the solution of the differential equation, the frequency 
of a tapered beam can be calculated if the solution for 
the corresponding uniform beam is known. 
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SYMBOLS 

6, angular displacement from equilibrium, 
y, z, rectangular coordinates in "beam. 
I-q , polar moment of inertia of section. 
J, torsion modulus of section. 
Gr, shear modulus of material. 
I , length of team . 
t , time. 

p, density of material. 

uu , angular frequency (<ju = 2tt£ , where f 
in cycles per second). 

U)} i angular frequency of a uniform team, 



TT / GJ 



21 y pip 

cp , torsion function. 
A, B , C x , C 2 , constants. 



Ill exponent in expression for variation of 

torsion modulus along span (n = 2X + l) 

X = (n - 1 )/2 . 

c , semi chord . 
/ \ 4/n 

5 = ? 



la) / ipP 1 OJ TT 



(i - r, t ) y j& i - Ti t wi 2 

t y j roots of equation (7). 
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m, number of mode. 

J x(H,Yx(t), Bessel functions of first and second kind, 

re spe c t i ve ly . 

Subscript s : 

o , at root. 

t , at t ip . 



ANALYSIS 



The differential equation for the motion of a section 
of a beam in torsional vibration (reference l) is 



p 1 * 



a 2 e 
at 2 



_a 

3z L 



GJ(z ) 



ae 



where GJ is the torsional rigidity of the section 
harmonic vibrations it is known that 



at 3 



- a) 2 e 



(i) 

For 



Then 



o r 



_d_ 

dz 



G J ( z ) 



d6 



p a) 2 L 9 = 0 



O + I M ii + x p p w2 " 6 = o 



dz' 



J dz dz 



JG 



(2) 



The functions I p and J (reference 4) for a prismatic 
"beam may "be expressed in the form 
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"D 



(x 2 + y 2 ) dx dy 



J = 



x 2 + y 2 + x — - y — J dx dy 
3y ox/ 



where cp , the torsion function, depends upon the cross 
section. This function is known for a number of shapes 
of cross section. A few examples from references 5 and 6 
are given : 



(a) Circle of radius 



r : 



^ = 0, J = I p = (rrr 2 ) 
(h) Ellipse with semiaxes a and b: 



a + "b a 2 + h 2 



(c) Bcctangle with sides 2a and 2h: 

(2n+l) ttx 



vtt/ n = o 



3 (~) n sinh 



2h 



, 3 (2n+l) Tra 

(2n+l ) cosh - 



s in 



( 2n+l ) ny 
2d 



2h 



J =3 — — k ao 

where k has the values given "by the following table 



a/b 


- ■ ■ ■ 

1.00 


1.50 


1.75 


2.00 


2.50 


3.00 


4 


6 


8 


10 


00 


k 


.424 


.589 


.642 


.668 


.748 


.790 


.844 


.89 6 


.922 


.940 


1 



The torsion modulus, J, of a section of a tapered 
"beam is very nearly the same as that of a prismatic beam 
having the same section. A comparison of the angular de- 
flections produced by a static torque in a conical shaft 
as calculated by the exact theory (reference 7) and by in- 
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tegrating the torsion modulus for a cylindrical shaft 
shows a discrepancy of the order of 1/6 tan 2 a, where a 
is the half angle of the cone. For a cone with tan a = 
0.10, the discrepancy is only 0.17 percent. In the fol- 
lowing analysis, this distinction in the torsion modulus 
has "been ignored. 

The "boundary conditions to he applied to the solution 
of equation (2) will depend upon the method of support of 
the "beam. For the case of most interest in this paper, a 
"beam built in at z = 0 and free at z * % $ the boundary 
condi t ions are : 



at z = 0 , 
at 



SB 

z = 1 ■ a 5 



6 = 0 

0 



A solution of equation (2) has "been found for the 
case of tapered "beams satisfying the following conditions: 

(a) All cross sections along the span are similar in 
shape . 

(o) The torsion modulus J can he expressed as a 
power of a linear function of position along the span. 
Then, from condition (a), the ratio tp/ J is a constant 

along the span and is equal to its value at the root sec- 
tion, I- 0 / J Q . From condition (h), 

- o 

J = J Q ( 1 - Bz ) n = J 0 ^ n 



The torsion modulus is proportional to the fourth power of 
the linear dimensions. Hence, in terms of the semichord, 

Tf = f = (- C ~) 4 
J o Vc o y 



4 /n 



At the tip, 



(3) 
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It may be noted that 



B = 



Figure 1 shews a number of typical plan forms corre- 
sponding to different values of n. The upper half of 
this figure can also he considered as a plot of semi chord 
against T| for different values of n. In this case T) 



increases from right to left 
these carves are noted: 



The following properties of 



at t| = 1, 
at % n 0, 



d c/c 0 n 
~d~~Tl ~ 4 

a c/ c 0 



d Tl 



== 0 when n > 4 
= oo when n < 4 
In terms of g and 3, equation (2) can he written 

For the case of a uniform "beam (dj/d| = 0), 
tion (2) has the well-known solution 

0 = 0! sin 3 uu I | + C 2 cos P U> I % 

The frequency equation for a cantilever beam is 

cos (3ui)l = 0 



e qua- 



P uu I - (2m ^ 1) 



TT 



where m = 1 , 2 , 3 , 
The lowest frequency ii 



OJ 1 



2pl 



Equation (2) can then be written 



d_6 + 1 d J d6 + 
d"p J d£ d| 



L TT ^ Q 



= o 



(4) 
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Now 

T] = 1 - B I | 



a - I ±1 - i_lJt 

6 " 1 - n t 

in a? V n V an 



d| aT) a| V V dTi 

After equation (4) has "been expressed in terms of T|, it 
"becomes 

2 2 2 

-n \ d 2 8 . ( -x n N n d6 , Aju_tt\ a - 0 



Put 



Then 



1 - T] t ai 1 2 0 



^ + n d_6 + e = (5 ) 

d£ 2 t d £ 



SOLUTION OF DIFFERENT IAL EQUATION 

A solution of equation (5 ) has "been found in terms of 
Bessel functions . Put 

8 - fx 'x(t) 

Then equation (5) "becomes 

d S J(t? + n_-_2* dJlli + A . ^---^U±li) J(t) 
. d£ 2 t N C ' 



1 
1* 



If X = ( n - l)/2, this equation reduces to Bessel r s equa- 
tion: 
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d 3 J(£) 1 dJ(£ ) A / A 2 



(l - p) J(t) - 0 



Therefore ( £) is a Bessel function of order A. 

The complete solution containing two arbitrary con- 
stants to be determined by the boundary conditions is 



9 = 



6 



7* 



i 



A^ (£) + 3Y A (£) 



X , non integral 



> (6) 



in tegral 



The function Y x (t) is a Bessel function of the second kind. 
The "boundary conditions in terms of £ are: 
at z = 0, t m L,« & = 0 



at 



d6 
d£ 



0 



After these conditions have been applied, the resulting 
equations can be simplified by means of the relations for 
Bessel functions: 



_ci_ 
dx 

_d 
dx 



x " J (x) = - x 
n j 



J . (x) 
n+ i 



x n J (x) 
n v ' 



xn , <x> 



The condition that the two equations for A and B have 
solutions other than A u 0, 3=0, is that the determi- 
nant of the coefficients of A and B vanish: 



WVo> 



J A + l^ t ^ 0 ) 



J -A^o 



"A 



~ J -A-l (Ti tU| 

^A(t 0 ) 
Y A + l^ t ^ 



1=0; X , nonintegral 



>- (?) 



0 ; X , integral 
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This condition determines an admissible set of values of 

t 0 , the roots of equation (?). These roots will he called 

tr- 

After the value or values of L that satisfy the 

det erminantal equation have "been found, the shapes of the 
deflection curves for the various modes can he directly 
obtained from equation (6). At the root section, £ = 

to = £r« 6 = °« 



0 = AJx (t r ) + BJ -x (t r ); 
0 = AJ X (l r ) + BY x (t r ); 



A , nonintegral 
X , integral 



Hence 
9 = 



A 

P 



-X b r 



] 



e . 3 
t* L 



*«>-|fe} MO] 



X , nonint egral 



integral 



V (8) 



where £ determines position along the span. 

Equation (7) is a transcendental equation of which 
the roots, £ r , determine the frequencies of the various 
modes. The equation cannot he solved explicitly "but may 
"be solved indirectly. The roots may he obtained to any 
desired accuracy by use of tables or by the graphical pro- 
cedure of plotting the value of the determinant against 
£ 0 . For each root there is a corresponding natural fre- 
quency of vibration given by 



?or the case of a pointed beam, equations (7) can be 
expressed in a simpler form. From the series expansion 
for Bessel functions, 



1 im 



J A+1 (\ V 



(yi t { fl) gfr+i) r(-x) m Q 



V 



-»o °-A-l 



lim 0 2(X+1) 



r (a +2) 



Hence, for pointed beams (Tl t = 0), the first of equations 
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( 7 ) do come s 



J% (t ) - - J X+l (1 Mo> J ft ) _ n 



When X is an integer, the sane result is obtained 

lim 

Tit 



— >o ^X+l^t^o' 



(\t 0 ) X+1 TT(T] t t 0 ) A+1 



= 0 



APPLICATION OF ANALYSIS 



The preceding analysis can "be immediately applied to 
teams the forms of which are mathematically specified by 
values of n and ffy . Most beams, however, will not ex- 
actly correspond to one of these cases; a method must 
therefore he found for determining the best values of n 
and for the actual "beam. An approximate method will 

he illustrated for a wing the shape of which is given hy a 
drawing. Other methods, such as the method of least 
squares, arc also available for determining these parame- 
ters. The thickness will he assumed to vary in the same 
way as the plan form. 

A sketch of a typical wing, together with a table of 
dimensions, is given in figure 2, It is required to rep- 
resent this shape by an equation of the form 

c = c 0 Tl 

The parameter T) varies linearly along the span from a 
value of T] t at the tip to 1 at the root. The procedure 
is to assume arbitrarily a value of \, then to make a 
logarithmic plot of chord or semichord against T). This 
process is repeated for different values of TL until the 
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value of T|x giving the best straight line is determined. 

The slope of this line is the value of n/4. The third and 
fourth columns of the table in figure 2 give values of T] 
for the assumed values for T| t of 0 and 0.1. Figure 3 

is a logarithmic plot of chord against T|. The points for 
T\ t = 0 fall nearly on a straight line. The slope of this 

line is 0.37, which gives n = 4 x 0.37 a 1.48. The values 
of c recalculated from the equation 

c ss 0.96 T] 0 ' 37 

are given in the last column of the table in figure 2. 

The frequency ratio, m/m x , determined from the roots 
of equation (7), is a function of only two variables, n 
and TL. Figure 4 shows this ratio as a function of Tit 
for several values of n. The actual frequency is given by 

<W X S 21 J I pP 

The frequency ratio for pointed beams (T] t = 0) is given 
in figure 5 for several values of n and for several modes. 
For the wing previously described (fig. 2), the frequency 
ratio is 1.77. The shape of the deflection curve is shown 
in figure 6 for several typical values of n and T\^. 

EXPERIMENTAL RESULTS 



Five duralumin beams (fig. 7) were constructed and 
their frequencies were measured to check the theoretical 
values of the frequency ratio. The central cross section 
of each beam was a rectangle 6 inches wide and l/4 inch 
thick. The beams were geometrically similar in cross sec- 
tion along the span and had straight taper, the correspond- 
ing values of rf t being 1.00, 0.75, 0.50, 0.25, 0.00. It 
is apparent from the boundary conditions that a freely sus- 
pended symmetrical beam vibrating in torsion with a node 
at the center will have the same frequencies as one built 
in at the center. The preceding analysis may therefore be 
applied to these beams. 



The experimental values of the frequency ratio arc 
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plotted together with the theoretical curve for n = 4 in 
figure 8. The greatest discrepancy "between experiment and 
theory is 2. percent, which is considered to he within the 
limit of accuracy determined "by the machining of the "beams 
and the measurements of the frequency. 



CONCLUSIONS 



The results of the tests of five "beams showed that 
the assumptions underlying the equation of torsional vi- 
bration of a tax>ered "beam 



at 2 



dz 



GJ ( z ) 



36 
3z 



are justifiable for all practical purposes. The 
ratio for even the most tapered case agreed with 
to within the experimental error. 



frequency 
the theory 



Langley Memorial Aeronautical Laboratory, 
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Figure 2.- Flan of a typical wing. 
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Figure 3.- Logarithmic plot of chord against r\ for a typical wing. 
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Figure 5.- The frequency ratio, UMBj, for pointed beams (r, t = 0) 
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Fig. 8 
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Figure 8.- Comparison of experimental values rf frequency ratio 
with theoretical curve for n« 4 # 



